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Abstract 

Two-phase random textures abound in a host of contexts, porous and composite media, ecolog- 
ical structures, biological media and astrophysical structures. Questions surrounding the spatial 
structure of such textures continue to pose many theoretical challenges. For example, can two- 
point correlation functions be identified that can be both manageably measured and yet reflect 
nontrivial higher-order structural information about the textures? We present a novel solution to 
this question by probing the information content of the widest class of different types of two-point 
functions examined to date using inverse "reconstruction" techniques. This enables us to show that 
a superior descriptor is the two-point cluster function 62(1*), which is sensitive to topological con- 
nectedness information. We demonstrate the utility of C*2(r) by accurately reconstructing textures 
drawn from materials science, cosmology and granular media, among other examples. Our work 
suggests an entirely new theoretical pathway to predict the bulk physical properties of random 
textures, and also has important ramifications for atomic and molecular systems. 

PACS numbers: 
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Two-phase random textures are ubiquitous in nature and synthetic situations. Examples 
include heterogeneous materials (e.g., composites, porous media and colloids) [l, 2], geologic 
media (e.g., rock formations) js], ecological structures (e.g., tree patterns in forests) 6], 
cosmological structures (e.g., galaxy distributions) j^, 5], and biological media (e.g., animal 
and plant tissue) 7| • Over a broad range of length scales, two-phase random textures exhibit 
a rich variety of structures with varying degrees of disorder and complex bulk properties 8^ 



10|, and questions concerning their quantitative characterizations continue to present many 



fundamental and practical challenges. 

It is well known that an infinite set of n-point correlation functions is generally required 
to completely statistically characterize such textures and their physical properties in the 
infinite- volume limit. A variety of different types of correlation functions arise in rigorous 
theories of structure/property relations One such basic quantity is the standard n-point 
correlation function S n (x.i, x 2 , . . . , x n ), which gives the probability of finding n points at 
positions xi, x 2 , . . . , x„ all in one of the phases Since the information contained in 

such an infinite set of functions is generally unattainable in practice, a natural starting point 
is to characterize the structure and bulk properties of random textures using lower-order 
versions. The two-point function S2, which is experimentally accessible via scattering of 
radiation 12], provides information about the distribution of pair separations. The three- 



point function S3 reveals information about how these pair separations involved in S2 are 
linked into triangles. The four-point function S4 controls the assembly of triangles into 
tetrahedra. A natural question is how much additional useful information do S3 and £4 
contain over and above S2? We will show that triangular and tetrahedral statistics do 
not significantly increase information content over and above pair statistics for textures 
possessing no long-range order. 

Therefore, an outstanding problem in condensed matter theory is to identify other two- 
point correlation functions that can be both manageably measured and yet reflect nontrivial 
higher-order structural information about the textures. The aim of this paper is to provide a 
novel solution to this problem using inverse techniques; specifically, "reconstruction" meth- 
ods. The purpose of a reconstruction (construction) technique is to reconstruct (construct) 
realizations of two-phase random textures that match limited structural information on the 
textures in the form of lower-order correlation functions 13], which are called the "target" 



functions. An effective reconstruction procedure enables one to generate accurate renditions 
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of random textures at will 
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20| and subsequent analysis can be performed on the recon 



HP 



struction to obtain desired macroscopic properties of the texture nondestructively 
Here we use this inverse methodology to determine the amount of structural information that 
is embodied in a set of targeted correlation function by quantifying the extent to which the 
original structure can be accurately reconstructed using those target functions. We quan- 
tify the accuracy of a reconstruction by measuring unconstrained (untargeted) correlation 
functions and comparing them to those of the original medium. 

We adapt the inverse reconstruction method of Yeong and Torquato jl^, Q to show that 
a superior two-point signature of random textures is the two-point cluster function C2(r) 



2l|, which is sensitive to topological connectedness information. We demonstrate that C^r) 



not only contains appreciably more information than 5*2, but more information than a variety 
of other "two-point" quantities, including the surface-surface correlation function F ss , the 
surface-void correlation function F sv , the pore-size function F, lineal-path function L, and 
the chord-length density function p (All the aforementioned correlation functions are 
defined precisely below.) Such a comprehensive study that investigates and compares all 
of the aforementioned statistical descriptors via inverse techniques has heretofore not been 
undertaken. Our results have practical implications for materials science, liquid- and solid- 
state physics, biological systems, cosmology, hydrology and ecology. 



I. DEFINITIONS OF THE TWO-POINT CORRELATION FUNCTIONS 

In the theory of random media, a variety of different types of two-point correlation func- 
tions naturally arise {if. Here we define seven different ones that will be employed in this 
paper. Consider d- dimensional two-phase textures in which phase i has volume fraction <pi 
(i = 1, 2) and is characterized by the indicator function 

*«<.) = ( 11 X£V " a) 

0, otherwise, 

where Vi is the region occupied by phase i (equal to 1 or 2). Note that here "phase" is used 
in a general sense in that it can refer to a solid, liquid or even void. The standard two-point 
correlation function is defined as 

^ ) (x 1 ,x 2 ) = <X«(x 1 )X«(x 2 )>, (2) 
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where angular brackets denote an ensemble average. This function is the probability of 
finding two points xi and X2 both in phase i. Henceforth, we will drop the superscript 
"2" and only consider the correlation functions for the phase of interest. For statistically 
homogeneous and isotropic microstructures, which is the focus of the rest of the paper, two- 
point correlation functions will only depend on the distance r = |x x — x 2 | between the points 
and hence S , 2(xi,x 2 ) = Sz{r). 

The surface- void F sv and surface-surface F ss correlation functions are respectively defined 

as 

F sv (r) = (M( Xl )X(x 2 )} , F ss (r) = (M( Xl ).M(x 2 )) , (3) 

where .M(x) = |VX(x)| is the two-phase interface indicator function. By associating a 
finite thickness with the interface, F sv and F ss can be interpreted, respectively, as the the 
probability of finding X! in the "dilated" interface region and x 2 in the void phase and the 
probability of finding both xi and x 2 in the "dilated" interface region but in the limit that 
the thickness tends to zero lj. 

The lineal-path function L(r) is the probability that an entire line of length r lies in 
the phase of interest, and thus contains a coarse level of connectedness information, albeit 
only along a lineal path The chord-length density function p(r) gives the probability 

associated with finding a "chord" of length r in the phase of interest and is directly pro- 
portional to the second derivative of L(r) 



. (Chords are the line segments between the 



intersections of an infinitely long line with the two-phase interface.) The pore-size function 
F(5) is related to the probability that a sphere of radius r can lie entirely in the phase of 
interest lj and therefore is the three-dimensional "spherical" version of the lineal measure 
L. 

The two-point cluster function C^r) gives the probability of finding two points separated 



by a distance r in the same cluster of the phase of interest 



211 ]. Note that a cluster of a 



phase is any topologically connected subset of that phase. The two-point cluster function 
can be measured experimentally using any appropriate three-dimensional imaging technique 
(e.g., tomography, confocal microscopy and MRI) [l|. The fact that C2 contains intrinsic 
three-dimensional topological information is to be contrasted with 5*2, which can obtained 
from a planar cross-section of the texture. In general, C2 is expected to embody a much 
greater level of three-dimensional connectedness information than either L or F, but the 
degree to which this is true has yet to be quantitatively demonstrated, which is one of the 



aims of this paper. 



II. INVERSE RECONSTRUCTION TECHNIQUE 

The stochastic optimization reconstruction algorithm for digitized media formulated by 
Yeong and Torquato 3] is ideally suited to carry out the aforementioned analysis because it 
can incorporate different t ype s o 



target statistical descriptors. This algorithm is both robust 
and simple to implement [15|-|20|. In this method, one starts with an initial realization of a 
random medium and a set of target correlation functions /*(R), /^(R), /n(R), • • •, which are 
obtained (i.e., measured) from the medium of interest. Here /"(R) is an n-point correlation 
function of type a, R = r 1; r 2 , . . . r„, and r, denotes the position vector of the ith point. The 
method proceeds to find a realization in which calculated correlation functions /*(R), /^(R), 
/^(R), . . . best match the corresponding target functions. This is achieved by minimizing 
an "energy" 

^EE[/"( R )-S( R )] 2 - ( 4 ) 

R a 

which is defined to be the sum of squared differences between the calculated and target 
functions, via a simulated annealing method in which a sequence of trial realizations is 
generated and accepted with the probability min{exp(— AE/T), 1}, where AE is the energy 
difference between the new and old realizations and T is a fictitious temperature. The 
initially "high" temperature is lowered according to a prescribed annealing schedule until 
the energy of the system approaches its ground-state value within a very small tolerance 
level. 



A. Degeneracy of Ground States Using S2 Alone 

It is instructive to consider the energy defined by (jl]) when only the standard two-point 
correlation function S2 is used. In this special instance, the energy is given by 

2 



E 



£ [S a (r) - S a (r)] , (5) 



where 5*2 and 5*2 are the two-point correlation functions of the target and reconstructed 
medium, respectively. We note that most previous reconstruction studies have only tried 
to match S 2 . However, it is now well established that ^(r) is not sufficient information to 



generally get an accurate rendition of the original microstructure |13Ml5l Il9l |20| In other 
words, the ground states when only S2 is incorporated in the energy [Eq. (jSJ)] are highly 
degenerate due to the nonuniqueness of the information content of this two-point function, 
which is clearly illustrated by the subsequent examples in the paper. The reader is also 
referred to the Supporting Text, which provides a rigorous explanation for the degeneracy of 
the ground states for digitized representations of textures. 



B. Insufficiency of Conventional Three-Point Information 

An obvious additional set of correlation functions that could be incorporated in the re- 
construction is the higher-order versions of S2, namely, S3, S4, etc. However, not only is the 
three-point correlation function S3 more difficult to compute, it is not at all clear that its in- 
corporation will result in appreciably better reconstructions because it only introduces local 
information about triangles when there is no long-range order, the most common occurrence. 
We can quantitatively verify the insufficiency of conventional triangular information by re- 



constructing a one-dimensional equilibrium distribution of equal-sized hard rods |24| using 
S3 (see the Supporting Text for technical details). Figure [T] compares this reconstruction to 
those involving £2 alone and a combination of S2 and C<i. It is clear that the reconstruction 
using 5*2 only significantly results erroneously in a highly clustered "rod" phase. Although 
incorporating S3 provides an improved reconstruction, it still contains large clusters and 
isolated rods that can be much smaller than actual rod size. On the other hand, the S2-C2 
hybrid reconstruction produces the most accurate rendition of the target medium. While 
this one-dimensional example is suggestive that C2 contains nontrivial structural information 
in excess to what is contained in S2, one must investigate this problem in higher dimensions, 
which presents algorithmic challenges, as we will describe in the next section. 



III. EFFICIENT ALGORITHMIC IMPLEMENTATION OF THE GENERAL 
PROBLEM 

Here we present a general methodology that enables one to efficiently incorporate a wide 
class of lower-order correlation functions in the Yeong-Torquato reconstruction procedure. 
The aforementioned probabilistic interpretations of the correlation functions enable us here 
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to develop a general sampling method for reconstruction of statistically homogeneous and 
isotropic digitized textures based on the "lattice-gas" formalism, which was introduced in 



Ref. 



20[ and has been generalized here. In the generalized formalism, pixels with different 



values (occupying the lattice sites) correspond to distinct local states and pixels with the 
same value are considered to be "molecules" of the same "gas" species 20|. The correlation 
functions of interest can be obtained by binning the separation distances between the selected 
pairs of molecules from particular species. 

In the case of S2, all "molecules" are of the same species. We denote the number of 
lattice-site separation distances of length r by Ng(r) and the number of molecule-pair sep- 
aration distances of length r by Np(r). Thus, the fraction of pair distances with both 
ends occupied by the phase of interest, i.e., the two-point correlation function, is given by 
S2O*) — Np(r)/Ns(r). To obtain C2, one needs to partition the "molecules" into different 
subsets Tj ("species") such that any two molecules of the same species are connected by 
a path composed of the same kind of molecules, i.e., molecules that form a cluster, which 
is identified using the "burning" algorithm 25|. The number of pair distances of length r 
between the "molecules" within the same subset Tj is denoted by N P (r). The two-point 
cluster function is then given by C2(r) = N P (r)/Ng(r). The calculation of F ss and F sv 
requires partitioning the "molecules" into two subsets: the surface set his containing only 
the "molecules" on the surfaces of the clusters and the volume set Ky containing the rest. 
In a digitized medium, the interface necessarily has a small but finite thickness determined 
by the pixel size. Thus, the surface-surface and surface-void correlation functions can be 
regarded as probabilities that are given by F ss = N ss (r)/Ng(r) and F sv = N sv (r)/Ng(r), 
respectively; where N ss (r) gives the number of distances between two surface molecules with 
length r and N sv is the counterpart for pairs with one molecule on the surface and the other 
inside the cluster. 

The lineal path function L can be obtained by computing the lengths of all digitized 
line segments (chords) composed of pixels of the phase of interest, and for each chord incre- 
menting the counters associated with the distances equal to and less than that chord length 
HQ. The chord-length density functions can then be easily obtained by binning the 



np 

tf 



chord lengths that are used to compute L The pore-size function F can be computed 
by finding the minimal separation distances of pixels within the phase of interest to those at 
the two-phase interfaces. The minimal distances are then binned to obtained a probability 
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distribution function, the complementary cumulative distribution function of which is F jsj. 

We have also devised methods to track clusters and surfaces that enable one to quickly 
compute the correlation functions of the new realization based on the old ones (see the 
Supporting Text for technical details) and thus make the Yeong-Torquato reconstruction 
procedure much more efficient than methods that directly re-sample the correlation functions 
for each trial realization. 

We have used this general procedure here to reconstruct a wide spectrum of random 
textures, including model microstructures, such as the cherry-pit model, equilibrium hard 
spheres, Debye random media and symmetric cell materials as weu as the digitized rep- 



resentations of sandstones |3|], metal-ceramic composites [26j, concrete microstructures [27]] , 
laser speckle patterns [2^] and galaxy distributions. Using the aforementioned largest set 
of correlation functions utilized to date, our analysis reveals that the best reconstructions 
always incorporate the two-point cluster function C^fV). In what follows, we will present 
specific results for only a subset of the correlation functions that we used, namely, various 
combinations of S2, F ss , F and C2, for the non-percolating phases of a concrete microstruc- 



ture [27|, a distribution of galaxies, and a three-dimensional hard-sphere packing. 



A. Concrete Microstructure 

The wide range of structural features in concrete, from nanometer-sized pores to 
centimeter-sized aggregates, makes it a wonderful example of a multi-scale microstructure 



281 ] . Figure [2]Ja) shows a binarized digitized image of a concrete sample cross-section. We 
have thresholded the original image so that the blue phase represents the stones and the 
lighter gray phase is the cement paste. The "stone" phase is characterized by a dense dis- 
persion of "particles" of various sizes: a nontrivial situation to reconstruct. Using S2 alone 
overestimates clustering in the system and indeed incorrectly yields a percolating "particle" 
phase. Thus, although S2 of the reconstruction matches the target one with very small 
error (see the figure in the Supporting Text), such information is insufficient to get a good 
reconstruction. Incorporating both S 2 an d surface- surface function F ss leads to a better 
rendition of the target system but the reconstruction still overestimates the degree of clus- 
tering. On the other hand, incorporating C2 yields an excellent reconstruction in that the 
"stone" phase clearly appears as a particle dispersion with a size distribution that closely 
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matches that of the target structure. As noted in the introduction, we can quantitatively 
test the accuracy of the reconstructions by measuring unconstrained correlation functions 
and comparing them to the corresponding quantities of the target system. Here we choose to 
compute the unconstrained lineal-path function L. Figure [2](e) reveals that the lineal-path 
function of the reconstruction that incorporates matches the target function L well and 
it clearly is appreciably more accurate than the other reconstructions. 



B. Galaxy Clusters 

Correlation functions have been used to understand the formation of galaxies and the 
large scale structure of the Universe j^, 5j. Such characterizations are becoming increas- 
ingly important with the advent of high-quality surveys in cosmology. We suggest that the 
reconstruction procedure can provide an important tool in cosmological and astrophysical 
applications, especially in helping to determine the lower-order correlation functions that 
reflect a priori information about nontrivial structural features, such as multi-scale cluster- 
ing and filamentary structures in the Universe. Figure [3]^a) shows a binarized image of a 
portion of the Abell 1689 galaxy cluster. We have chosen the binarizing threshold such that 
the "galaxy" phase (bright spots in Fig. |3k) exhibits different sized clusters. As one can see, 
the "primary" (single largest) cluster in the target system has been reproduced by all the 
reconstructions. This is because the volume fraction of the "galaxy" phase is relatively small 
while the target S2 has a relatively long tail (see the figure in the Supporting Text), which 
requires the clustering of a majority of "galaxy" phase. However, the "secondary" (smaller 
compact) clusters are significantly different for different reconstructions. The reconstruction 
using S2 only produces large elongated secondary clusters. Incorporating F ss enables one to 
obtain a better rendition, however it still contains a few elongated secondary clusters. The 
incorporation of C2 again provides the most accurate reconstruction. This is visually verified 
by examination of both the sizes and shapes of the primary and secondary clusters, and is 
quantitatively confirmed by the comparison of the unconstrained L of the reconstructed and 
target systems, as shown in Fig. [3te). 
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C. Sphere Packing 



As an application of our methodology to three dimensions, we have reconstructed a 
digitized realization of an equilibrium distribution of equal-sized hard spheres, as shown in 
Fig. Ufa). This packing is generated using the standard Metropolis Monte Carlo technique 
: or a canonical ensemble of hard spheres in a cubical box under periodic boundary conditions 
l| . A visual comparison of the hybrid reconstruction involving the two-point cluster function 
reveals that it accurately yields a dispersion of well-defined spherical inclusions of the same 
size, in contrast to the S2 reconstruction, which again grossly overestimates clustering of 
the "sphere" phase. In contrast to the previous examples, here we incorporate the pore-size 
function F (not the surface correlation functions) in one of the reconstructions. Although 
the reconstruction incorporating F provides improvement over the rendition on the SValone 
reconstruction, it is still inferior to the S2-C2 reconstruction in reproducing both the size 
and shape of the "sphere" phase. (The target and reconstructed correlation functions are 
shown in the Supporting Text.) The accuracy of the S2-C2 hybrid reconstruction can also be 
seen by comparing the unconstrained lineal-path function L of the target system to those of 
the reconstructed media [see Fig. BJe)]. 



IV. DISCUSSION 



In summary, while it was known that the information content of the standard two-point 
function S2 of a random texture is far from complete, we have demonstrated here that 
the next higher-order version S3 generally does not contain appreciably greater informa- 
tion. The fact that this natural extension to incorporate higher-order S n , which has been 
pursued in the last century in statistical mechanics, is not a fruitful path motivated us to 
inquire whether there exist sensitive two-point statistical descriptors that embody nontriv- 
ial structural information. We probed the information content of seven different types of 
two-point functions using inverse reconstruction methods. For all of the examples studied 
here, reconstructions that include the two-point cluster function C2 were always found to 
be significantly more accurate than those involving any of the combinations of pairs of the 
other functions. 

More precisely, the incorporation of C2 significantly reduces the number of compatible 
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microstructures as compared to the compatible microstructures consistent with the same 
three-point function S3, which is schematically indicated in Fig. [5j The two-point cluster 
function is an especially sensitive structural signature when clustering and phase connected- 
ness are present, the most difficult situations to treat. This can be seen from the reconstruc- 
tion of the hard- rod system that incorporated C 2 (see Fig. [I]). More importantly, we showed 
the utility of 62(1") in higher dimensions by accurately reconstructing galaxy distributions, 
concrete microstructures and dense hard-sphere packings, among other examples. 

Why is C 2 (r) a superior two-point structural signature? To answer this question, it is 
useful to first compare it to S 2 (r). The latter, unlike C 2 (r), does not distinguish between 
events in which the end points of the line segment of length r fall in the same cluster 



2l| 



of a particular phase and those that do not involve the same cluster of that phase 
More precisely, C 2 is the "connectedness" contribution to the standard two-point correlation 
function, i.e., 

S 2 (r) = C 2 (r) + D 2 (r), (6) 

where D 2 measures the probability that the end points of a line segment of length r fall in 
different clusters of the phase of interest. Therefore, whereas S 2 (r) is insensitive to clustering 
and percolation, C 2 (r) becomes a progressively longer-ranged function as clusters grow in 
size such that its volume integral diverges at the percolation threshold . By contrast, 

the quantities L, p, F, F ss and F sv are insensitive to crossing the percolation threshold. 
Indeed, for particle systems, one can show that C 2 is a functional of the infinite set of 



"n-particle" connectedness functions 21]. 



Thus, although C 2 is a "two-point" quantity, it actually embodies higher-order structural 
information in a way that makes it a highly sensitive statistical descriptor over and above S 2 . 
This ability to "leap frog" past the usual approach of incorporating additional information 
via higher-order versions of S 2 has important ramifications for new structure/property of 
random textures. Specifically, our work suggests that new theories should be developed 
that relate the transport, mechanical, chemical and optical properties of random textures 
to functionals that incorporate C 2 . It is clear that such theories would be highly predictive, 
since it is well established that the presence of clusters in textures can dramatically alter 
their macroscopic physical properties [l|. 

It should not go unnoticed that our work also has important implications for atomic and 
molecular systems where the analogous standard correlation functions that arise are the 
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two-body correlation function g 2 , three-body correlation function g%, etc. 29M32| The three 



body function g 3 , for example, has been the focus of great attention, e.g., integral equations 



and approximations have been devised for g 3 [30b32|] . The spatial structure of disordered 
atomic or molecular systems may be regarded to be special cases of random textures. For 
example, the particles that comprise simple atomic systems are fully specified by their center- 
of-mass coordinates. These point distributions can be decorated in an infinite number of 
ways to produce random textures in the sense that we have defined them in this paper. For 
example, one could circumscribe the points by spheres of size dictated by the physics of 
the problem (e.g., electron hopping distance). For such systems, our work suggests that the 
pair-connectedness function P 2 (i.e., the connectedness contribution to g 2 [l|) should contain 
far greater information than g 3 beyond that contained in g 2 . 
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FIG. 1: (a) Target: an equilibrium hard-rod system in which the equal-sized rods cover 50% space. 
Each rod in the system is 10 pixels in length, (b) 52-alone reconstruction, (c) S3 reconstruction, 
(d) S2-C2 hybrid reconstruction. For visualization purposes, the one-dimensional rod systems are 
artificially extended in the vertical direction. 
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FIG. 2: (a) Target system: a binarized image of a cross-section of concrete [27[. The linear 
size of the digitized texture is Nl = 170 pixels, (b) i>2-alone reconstruction, (c) S<z-F ss hybrid 
reconstruction, (d) S2-C2 hybrid reconstruction. All the reconstructions are associated with a final 
energy (error) E ~ 10 -8 . (e) The unconstrained lineal-path function L of the reconstructions and 
the target image. Pixel size supplies the unit for the distance r. 
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FIG. 3: (a) Target system: a portion of the Coma Cluster. The linear size of the digitized texture 
is Nl = 238 pixels, (b) SValone reconstruction, (c) S2-F ss hybrid reconstruction, (d) S2-C2 
hybrid reconstruction. All the reconstructions are associated with a final energy (error) E ~ 10~ 8 . 
(e) The unconstrained lineal-path function L of the reconstructed and target systems. Pixel size 
supplies the unit for the distance r. 
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FIG. 4: (a) Target system: a digitized realization of a hard-sphere packing in which the spheres 
occupy 44.6% of space. The linear size of the digitized texture is Nl = 100 pixels, (b) S^-alone 
reconstruction, (c) S2-F reconstruction, (d) S2-C2 hybrid reconstruction. All the reconstructions 
are associated with a final energy (error) E ~ 10 -11 . (e) The (unconstrained) lineal-path function 
sampled from target and reconstructed realizations. Pixel size supplies the unit for the distance r. 
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FIG. 5: The set of all micro-structures associated with a particular S2 is schematically shown as 
the region enclosed by the solid contour in both panels. The shaded region in the left panel shows 
the set of all microstructures associated with the same S2 and S3. The shaded more restrictive 
region in the right panel shows the set of all microstructures associated with the same £2 and C2. 
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